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1.  INTRODUCTION:  NEW  SHIP  STRUCTURE  AND  MODEL 

The  proposed  new  surface  ship  hull  concept  consists  of  a  double  skin  that  wraps 
around  the  bottom,  sides,  and  main  deck.  The  two  skins  are  connected  by  plates  normal 
to  the  surfaces  forming  a  cellular  structure  similar  to  a  cardboard  box,  as  shown  in  Figure 
1.  The  advantages  of  this  type  of  construction  have  been  described  by  Okamoto,  et.  al. 
[1],  and  Beach  [2]. 

The  object  of  this  report  is  to  develop  formulas  for  estimating  the  principal  stresses 
in  such  a  structure.  In  order  to  be  able  to  obtain  simple  analytical  formulas,  we  model  a 
ship  hull  as  a  circular  cylindrical  orthotropic  shell  surrounding  an  elastic  core.  Under  a 
pure  end  bending  moment  M  and  a  lateral  pressure  loading  q  which  does  not  vary  along 
the  length,  the  cylinder  bends  into  a  curved  tube  of  oval  cross  section,  as  shown  in  Figure 
2. 


2.    SHELL  EQUATIONS  AND  SOLUTION 

The  theory  governing  the  deformation  of  thin  shells  is  well  developed  and  available  in 
many  forms.  Here  we  use  the  semi-momentless  shell  theory  of  Axelrad  [3].  **  When  spe- 
cialized to  the  linear  St.  Venant  problem  for  circular  cylindrical  shells,  Axelrad's  equations 
(2.122)  -  (2.123)  reduce  to 

dr)4    +   dif  D2  drj*  K  ] 

V  +  V=°  (2) 

Here  the  coordinate  r\  denotes  the  circumferential  angle  and  R  denotes  the  radius  of  the 
undeformed  cylinder,  as  shown  in  Figure  2.  The  principal  values  of  the  change  in  curvature 
tensor  are  denoted  by  /vj (77)  and  k2(tj) ;  their  geometrical  definitions  are  illustrated  in 
Figure  3.  The  longitudinal  membrane  force  ^(77)  (force  per  unit  length  of  the  shell 
midsurface)  is  related  to  the  longitudinal  strain  e^r?)  by  the  constitutive  law 

Ti  =  Ehe2  (3) 

Here  h  is  the  thickness  of  the  shell,  and  E  is  Young's  modulus.  The  longitudinal  bending 
moment  M\(r))  and  transverse  bending  moment  M2(rj)  (moments  per  unit  length  of  the 
shell  midsurface)  are  related  to  the  curvature  changes  «i(ry)  and  k2(tj)  by  the  constitutive 
laws 

Mi  =  D^Ki  +  uk2)  >     M2=D2(k2  +  vki)  (4) 


** 


Our  equations  may  also  be  obtained  from  other  shell  theories,  such  as  that  of  Sim- 
monds  [4]. 


Here  Di    and  D2   are  bending  stiffnesses,  and   v  is  Poisson's  ratio.    The  longitudinal 
curvature  change  «i(ry)  and  transverse  membrane  force  Ti{r\)  are  given  in  terms  of  the 

above  variables  by 

1     d2Tx  1  d2M2  ,  x 

Kl  =  -EhR~W'     T2=Rq+R~d^~  (5) 

The  net  bending  moment  M  acting  on  any  cross  section  of  the  tube  may  be  calculated 

from  the  formula 

r2* 

M  =  R2  Tjcosry  dr)  (6) 

Jo 

The  strain  measures  are  related  to  the  displacement  components  by  the  equations 

1  d2w 

1    (d2w        \ 

K*=-&W+W)  (8) 

_1_  /  d2w       Zdv      ldu\ 
T~     R*\dtdri      4d£  +  4^,)  (9) 

Idu 
«.-5S£  '       (10) 

1  ( du      dv\  .     . 

7-5U  +  3f)"°  (12) 

Here  R^  denotes  the  distance  measured  along  the  axis  of  the  undeformed  tube,  as  shown  in 
Figure  2,  and  (u,  t>,  w)  denote  displacement  components  in  the  (f  ,r? ,  radial)  directions. 
We  suppose  that  the  pressure  loading  q(rj)  acting  normal  to  the  surface  of  the  shell 
can  be  expressed  as 

en 

q{ri)  =  -K[w(r))  +  w{r)  +  tt)]  -  —(1  +  cos  2  77)  (13) 


at  rj  =  0  is  the  longitudinal  curvature  of  the  cylinder  axis  due  to  the  beam-like  bending. 
The  terms  in  (15),  (17),  (19),  (20)  which  are  underlined  are  secondary  quantities  arising 
from  ovalization  of  the  cross  section.  Proper  design  requires  that  these  quantities  be  made 
relatively  small  by  having  a  large  enough  foundation  modulus  K . 


3.   APPLICATION  TO  DOUBLE  HULL 

As  an  example  of  the  application  of  these  formulas,  let  us  consider  the  double  hull 
sketched  in  Figure  1.  An  average  hull  section  contains  3  plates  each  of  thickness  t  and 
width  b ,  as  shown  in  Figure  5,  so  the  average  hull  thickness  is 

h  =  3t  (21) 

The  bending  stiffnesses  corresponding  to  this  double  hull  geometry  are 

7Eb2t  „  Eb2t 

D>  =  w^r  d'  =  w^t-  (22) 

Of  interest  for  design  purposes  are  the  maximum  stresses  arising  in  the  double  hull.  The 
magnitudes  of  the  membrane  stresses  corresponding  to  the  membrane  forces  (14)  -  (15) 

are  less  than 

M 
°Ml  =  MR2  (23) 

8R2  ub2M  6R2 

°M2~    4*    +12*(l-,/»)*iy  + 121(1  +  ^)  (24) 

The  magnitudes  of  the  bending  stresses  corresponding  to  the  bending  forces  obtained  from 

(4)  and  (16)  -  (17)  are  less  than 

bM  u6R* 

aBl       7tt(1  -  „*)iJ*3  +  76t(1  +  ^gl)  t25) 

6R3  vbM 

°m         66/(1  +2£Kl)  +  67r(l-  v2)tR3  (26) 

FVom  the  architextural  drawings  of  the  double  hull  cross  section  sketched  in  Figure 
1,  we  have  computed  the  average  values  displayed  in  the  Table  of  the  various  geometrical 
and  material  parameters  in  our  formulas.  The  value  of  M  is  the  maximum  value  of  the 


design  hogging  bending  moment.  I  is  the  moment  of  inertia  of  the  hull  cross  section  about 
its  centroid.  If  we  choose  the  cylindrical  cross  section  to  have  the  same  moment  of  inertia 
as  the  hull,  this  determines  the  radius  of  the  cylinder  to  be  R  =  25.75  ft.  We  suppose  that 
the  foundation  stiffness  K  is  the  same  as  the  bulkhead  stiffness.  From  experiments  on 
pyramidal  truss  cores  typical  of  those  which  are  proposed  for  the  bulkhead  structure,  we 
have  determined  the  average  value  of  K  shown  in  the  Table. 

Now  let  us  put  some  of  the  numbers  shown  in  the  Table  into  the  formulas  (23)  -  (26): 

aMi  =  21.6  (27) 


1  "*"     9D2 


afii  =  1-15    +    ;  ,  2kr<  (29) 

1  +     9D2 


12 A 

9D2 


^=n    2KR<         +        A  (30) 


Here  all  stresses  are  in  units  of  ksi.  Note  that  if  there  were  no  elastic  foundation  (K  =  0) , 
the  secondary  stresses  (those  terms  underlined  in  (28)  -(30))  would  be  a  sizeable  fraction  of 
the  primary  stresses.  But  gffi  =  3733 ,  from  the  values  in  the  table,  so  these  secondary 
stresses  are  rendered  negligible  by  the  bulkhead  stiffness. 


4.    CONCLUSION  AND  FUTURE  WORK 

Our  analytical  formulas  indicate  that  the  stiffness  of  the  proposed  bulkheads  could 
be  reduced  by  a  factor  of  10  without  incurring  significant  secondary  stresses  in  the  double 
hull.  A  less  bulky  bulkhead  design  would  have  obvious  cost  and  weight  benefits. 

This  analysis  models  only  the  most  fundamental  aspects  of  the  complex  state  of  stress 
which  could  actually  exist  in  the  hull  of  a  ship  at  sea.  In  future  work  the  model  should  be 
refined  to  make  its  predictions  correspond  more  closely  to  reality.  Features  which  should 
be  included  are: 

1.  A  more  rectangular  cross  section  corresponding  to  the  actual  hull  shape.  This  could 
perhaps  be  included  within  the  framework  of  the  present  analysis  by  conformal  map- 
ping techniques. 

2.  A  more  realistic  applied  loading,  including  torsion  and  internal  pressures. 

3.  Discrete  treatment  of  the  bulkheads  and  other  stiffeners.  Also,  allowance  for  nonuni- 
form stiffness  properties. 

4.  Geometric  and  material  nonlinearities.  This  is  essential  for  an  ultimate  strength  anal- 
ysis. 
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Figure  1:    Double  hull  designed  by  David  Taylor  Research  Center. 
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Figure  2:    Idealized  Model  of  Double  Hull 
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Figure  3:    Geometrical  Definitions  of  Curvature  Changes 
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q  =  -—  (1  +  cos2t]) 


Figure  4:    Lateral  Pressure  Loading 
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Figure  5:    Stresses  Acting  in  Double  Hull  Wall 
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Table  1:    Typical  Parameter  Values 


M 

665,000  kip-ft 

6 

0.0624  kip/ft3 

E      ^^H 

29,500  ksi 

V 

0.3 

t 

0.41  in 

b 

35  in 

I=7iR3h 

5500  ft4 

K 

15  kip/in3 
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